Optimal Quantum States for Image Sensing in Loss 
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We consider a general image sensing framework that includes many quantum sensing problems 
by appropriate choice of image set, prior probabilities, and cost function. For any such problem, in 
the presence of loss and a signal energy constraint, we show that a pure input state of light with the 
signal modes in a mixture of number states minimizes the cost among all ancilla-assisted parallel 
strategies. Lossy binary phase discrimination with a peak photon number constraint and general 
lossless image sensing are considered as examples. 

PACS numbers: 42.50.Ex, 42.50.Dv, 06.20.-f, 03.67.Hk 
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The use of nonclassical and entangled states of light, 
i.e., states other than the easily generated coherent states 
and their classical mixtures [T], for applications such 
as sub-shot-noise imaging [51 and imaging with sub- 
Rayleigh resolution ^ has received much attention. In 
the areas of sensing and metrology [3] , there have been re- 
cent theoretical studies of quantum-enhanced target de- 
tection [S] , reading of a digital memory [SI [7] , and of op- 
tical phase estimation [S] with nonclassical states. Given 
the interest in applications of quantum states of light for 
sensing, it is important to theoretically establish what 
state(s) accomplish a sensing task using minimum en- 
ergy, assuming the most general measurements and post- 
processing. This would place a limit on the enhancements 
obtainable from nonclassical states using experimentally 
realizable measurements. The ubiquitous linear loss is 
known to be a bottleneck for harnessing quantum advan- 
tage in many communication and metrology applications 
[SI [3] . Although the problems of [SHT] naturally include 
various degrees of loss, few general results including its ef- 
fects are available. In this Letter, we first set up a general 
framework for image sensing in the presence of loss that 
subsumes many of the above problems. We then identify 
a class of input states that contains an optimal, i.e., cost- 
minimizing, state for any problem fitting the framework, 
and under any form of signal energy constraint. 

General Image Sensing Framework: Suppose an im- 
age is drawn, unknown to the receiver, from a set X = 
{/i, . . . ,/m} of M images according to the probability 
distribution {tti, . . . ,ttm}- We model each image as a 
pixelated (transmissive or reflective) optical mask with 
uniform transmissivity/reflectivity and phase shift within 
each pixel. For P the number of pixels in each image, the 
p-th pixel (p g {!,..., P}) of image /„ is modeled as a 
beam splitter effecting the mode transformation 
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Here rjm' is the transmittance (or reflectance if reflec- 
tive probing is used) of the p-th pixel in Im, and 9m is 
the phase shift imparted to the input (or "signal") field 
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FIG. 1. Schematic of procedure for sensing of an unknown 
image from a set X with pixels described by (r;™ , 0m ), p = 
1,...,P, via ([T|. The source generates signal modes {af} 
for probing the image and idler modes which are retained 
losslessly. An optimal measurement for a; G A" is made on the 
J return modes {bf } and J' idler modes jointly. 



modes probing the p-th pixel of /,„ alone. 

For probing the unknown image, we consider quantum 
states of J = X]n=i -^^^^ signal modes that may be entan- 
gled to ancilla (or "idler") modes that are not sent out 
to interrogate the image but held losslessly (see Fig. 1). 
Here, J^^' is the number of modes interrogating pixel 
p, which, for example, could be successive time modes. 
Eq. ^ includes the annihilation operator dj of the j- 
th signal field mode probing pixel p and the annihilation 
operator ej of the j-th input environment mode at p. 
The input environment modes are taken to be in the vac- 
uum state - the assumption of no thermal noise in the 
environment is realistic at optical frequencies. We fur- 
ther assume that the output mode corresponding to 6^ 

in (fTl) , but not that corresponding to f^^' , is available for 
making quantum measurements. This is almost always 
the case in practice as the environment input and output 
modes are not easily accessible to the user in a standoff 
imaging scenario, or if the light source and receiver are 
at different spatial locations. Additional loss during state 
propagation may be included as a multiplicative factor in 
the {Tyi^)}. 

We first consider pure input states - we will return 
to the mixed input state case later. An arbitrary pure 



quantum state \tp)is of the signal and idler modes may 
be written in the form 
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Here n = {n\ , . . . , ,,jf,-^, . . . , ,,^ 

dimensional vector whose component n^ 
photon number in the j-th mode interrogating the p-th 
pixel, and |n)5 are Fock states of the signal modes. We 
do not restrict the number J' of the idler modes, nor the 
form of the idler states \(f)n)i, only requiring that they be 
normalized. Allowing an arbitrary input state (pi) corre- 
sponds to the most general ancilla- assisted parallel strat- 
egy (see H] for discussion on parallel strategies). Little 
is known about non-parallel strategies (but see [TU]) that 
may include adaptive selection of inputs which is known 
to assist some channel discrimination problems [11] . Irre- 
spective of the form of the (l^n)/}: the probability mass 
function (pmf ) of the photon number in the signal modes 
ispn = |cnP, which determines quantities of interest such 
as the mean total signal energy:- 

\p=l i=l / n 
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In practice, the mean total signal energy may be upper 
bounded by a given number Ns- 

Once an input state \^)is is chosen and the signal 
modes are sent to probe the image, the return-|-idler 
states constitute an ensemble £ = {{TTm, Pm)}m=i^ where 
Pm = id/ ^^m(|V')/s(V'l) is the density operator on the 
return-|-idler Hilbert space at the output of the quantum 
channel idj iSilCm resulting from the interaction of the sig- 
nal modes with /^ via M and the identity map on the 
idler modes. Depending on the imaging task, we attempt 
to extract a parameter lying in an observation space X 
by making a quantum measurement that is represented 
by a POVM [T^] with outcomes x ^ X and correspond- 
ing operators {Ex}xex- The task also specifies a cost 
function C{m,x), and we are interested in the minimum 
average cost C 
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where the minimization is over all POVMs {Ex}x£X- 
Therefore, adaptive measurements are included in our 
model. Note that the input state and image param- 
eters enter into the cost via the ensemble £ while the 
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FIG. 2. Image sensing problems in linear loss. Left: Quantum 
reading - Digital reader composed of a transmitter T and re- 
ceiver i?. Right: Interferometer for discrimination/estimation 
of the relative phase shift S with two-mode signal-only source 
S and detector D. 



imaging task determines X and the cost function. Thus, 
choosing X — {1, . . . , M} and C{m, x) ~ 1 — 5„i^x makes 
C[£] equal the minimum probability of error (MPE) in 
discriminating the M images. For the same cost func- 
tion, Af = 2, P = 1, and 9m = corresponds to the 
quantum reading and target detection problems of [5]- 
T. For P = 1 and rf^^ = t], choosing X = [0, 27r) 
and C{m,x) — [x — 9m ]'^ corresponds to minimum- 
mean-square-error (MMSE) phase discrimination in the 
presence of loss. As M — >■ oo, we approach MMSE 
phase estimation in loss. The usual interferometric setup 
for phase estimation [U [5] is recovered using P — 2, 
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and the cost 



m={S,(^)G [0,27r)2, .Tn - 
function C{m,x) = [x — 6]"^. Here 5 is the relative phase 
shift of interest and (j) is the undesired common phase 
shift in both arms of the interferometer (see Fig. 2). 

In the state ([2]), if It/fn)/ = |</>)/ for all n and some 
idler state | (/>)/, the state factorizes so that the idler 
is effectively absent. The other extreme is the case of 
(0n|0n')/ = Sa,n' SO that the density operator on the sig- 
nal modes is diagonal in the multimode Fock basis. Such 
states are called Number-Diagonal Signal (NDS) States 
in [7j , well-known examples being the two-mode squeezed 
vacuum of [3 [6] and the NOON state [13]. In [7], the 
error probability and other quantities of interest were 
computed for NDS inputs in the M = 2 case. Our main 
result is that the NDS states are optimally 'matched' to 
the general imaging problem described above. 

Theorem 1 (NDS State Lower Bound). Let I be 
a set of M P-pixel images described via transforma- 
tions of the form (fl]) with prior probabilities {TTm}m=i- 
For any imaging task with cost function C{m,x), the 
minimum cost C achieved by the input state |?/')/s = 
^j^ Cn|0n)/|n)5 is lowcr bounded by that achieved by 
a corresponding NDS state \^)is — J2n'^n\^n)i\'>^)sj 
where {|<&n)/} is any orthonormal set, with the same sig- 
nal photon pmf. 

General Results in Quantum Decision Theory: The 
proof of Theorem 1 requires two simple but general re- 
sults in quantum decision theory. To state them, we de- 
fine the notion of mixture of ensembles:- For each value 



of an arbitrary index I with associated probability A; , let 

£i = {(I'm , /Om )}^=i be an M-ary ensemble of states in 
a Hilbert space T-l. Then the M-tuple of pairs 
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is also an ensemble, called the mixture of the ensembles 
{£i} - the {£[} are sub-ensem,bles of £. A mixture of 
ensembles can arise from a two-step procedure in which 
I is chosen with probability A;, following which the state 
pirl is prepared with probability TTm ■ In our proof of 
Theorem 1, sub-ensembles arise as the conditional output 
states of a measurement with outcomes {/} on a given 
ensemble. We have the following basic result: 

Lemma 1 (Concavity of C under mixing of ensembles). 
Consider a sensing task with cost function C{m,x). For 
M-ary ensembles {£i} indexed by I, and probability dis- 
tribution {A;}, 
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The notion of orthogonal ensembles provides a sufficient 
condition for equality in (21). The support of an en- 
semble £ — {(tTto, /Om)}m=i is defined to be suppf := 
^^^j^ ranp„j, where ranpm is the range of p„i. Then 
ensembles £ and T are said to be orthogonal if the sup- 
port spaces suppi? and supp J^ are orthogonal. 

Lemma 2. If {£i\ are pairwise orthogonal ensembles on 
H, we have 
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J2XA =J2\iC[£i]. 
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The proofs of Lemma 1 and 2 are given in the Appendix 
along with their physical interpretation. 

Proof of Theorem 1: Suppose the general state (|2| is 
used as input. To calculate the output states {/Om}m=i, 
we may use the Schrodinger-picture form of (fll) to get 
the purification 

l'^-) -Hill c„A(-')|0„),|n - \)s\ \\)e (7) 
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where \\) e is a Fock state of the environment modes and 
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The output state Pm is then given by (note that \ipm)is 
are non- normalized states): 
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In (TM), 1 is the (random and unknown) pattern of the 
number of photons leaked from the signal modes into the 
environment modes during interrogation of the image. 
The probability that the leaked photon pattern is 1 is 
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so that the conditional probability of hypothesis m given 
lis 



TT ^ . 

Ai 

Thus, the ensemble 

£ = {{T^m, Pm)}m=l = X, ^1^1 

1 
for the sub-ensembles {£\} given by 
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According to Lemma 1, the mixture (12) satisfies 
C[£]>Y^ XiC [£i] . 
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Consider the RHS of (14 1. For each 1, £i is a pure-state 



ensemble, so C[£{\ is a function of just the pairwise inner 

products (V'£Vm')- = '^''' 

m-- 
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The crucial point is that, owing to the form (IT]) of the 
beam splitter transformation, G,^ „j, is independent of 

the choice of the {|</>n)/}- From ( lOpT ), so are Ai and 

(1) — — 

TTm . Thus, the hypothetical measurement scenario in 

which one has knowledge of 1 (or alternatively, one is 

allowed to make a photon number measurement on all 

the output environment modes) and whose C is given by 



the RHS of (14 1, has the same C for any choice of the 

{|0n)/}. 

Finally, we consider the NDS input state |^)/5 = 
J2n'^^\^n)i\^)s corresponding to ^ satisfying 

($„|$n')/='5n,n'. (16) 



It is readily verified using ^ and ( 16 1 that 

(V'^^l£''>/s = V(V'^^lV'^',)/s (17) 

so that the {£1} are pairwise orthogonal ensembles. 
Therefore, by Lemma 2, the NDS input |^)/s attains 
^[ AiC [£i\, and does so with the same signal photon pmf 
as |V')/5- _ 

Discussion and Implications: The C of the NDS state 
\'^)is of Theorem 1 is a function of only the J-mode 
photon pmf {pn}- Thus, for a given J, the search for an 
optimal input state for a given imaging task may be con- 
fined to the set of {pn} satisfying given constraints, e.g., 
an average/peak signal energy constraint or a mode-by- 
mode signal energy constraint. As illustrated below, the 
problem of finding the optimal quantum state reduces to 
the classical problem of finding an optimal probability 
distribution. To see that mixed input states pis do not 
help, we first purify p/5 using added idler modes. As 
(pi) contains no restriction on the idlers, the NDS state 
corresponding to the purification - which has the same 
{Pn} as Pis - has C not larger than that of the purifica- 
tion (which in turn beats pis)- Note also that the perfor- 
mance achieved by an arbitrary state of signal energy Ns 
can be achieved by an NDS state of total (signal-f-idler) 
energy not larger than 2Ns by choosing \(f>n)i — |n)/, 
the Fock state of the idler modes. 

Theorem 1 strongly suggests that ancilla-assisted par- 
allel strategies for image sensing are superior to signal- 
only parallel strategies. This is known for discrimina- 
tion between some pairs of channels [15 and is also true 
for our phase discrimination example below. We conjec- 
ture that they are strictly better whenever nonzero loss 
is present (77m < 1) because for a signal-only input state, 
the {£1} are not orthogonal and arc unlikely to achieve 
the lower bound of (14 1. Such ancilla-assisted schemes 



appear to be unexplored for some problems of interest - 
e.g., studies of the optimal state for phase estimation in 
loss have hitherto been confined to two-mode signal-only 
states [8]. At the same time, Theorem 1 implies that the 
best possible performance can be obtained without ancil- 
lary modes if the value of 1 is known. This result should 
place interesting limitations on the quantum advantage 
obtainable in any sensing problem. 

Binary Phase Discrimination: As an application of 
Theorem 1, we obtain the single-pass {J — 1) state that 
discriminates between a and tt phase shift with mini- 
mum error probability among states with a peak signal 
photon constraint of A'poak = 2 in the presence of loss. 
In the terminology of our framework, we have M = 2, 



P 



1, vi'' 



V^' 



ri <\, and 6\ 



(1) 



0, 



3(1) 



TT. We 



also assume vri = 7r2 = 1/2. An arbitrary state \ip)is 
satisfying these constraints is 



|V')/s = Vp^|0o>/|O)s + ^/P~i\<t>i)i\l)s + 



P~2\4>2)l\2)s, 

(18) 



where phase factors have been absorbed into the normal- 
ized kets {\4>n)i}- In terms of the density operators pi 
and p2 defined earlier, the minimum error probability is 
given by the Helstrom formula [H] 
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where || • ||i is the trace norm. According to Theorem 1, 
we may confine our search for optimal states to the NDS 
class for which {If/in)/} are orthonormal. For such states, 
the minimum error probability is given in closed form by 
Eq. (39) of [7]: 
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P, = 1/2 - ypi ^(pory + P2nr^' + (2P2^(1 - vY) 

Since pq +pi +p2 = 1, we may consider po and pi as inde- 
pendent variables taking values in the triangle T whose 
vertices have the [po,pi) values (0,0), (1,0), and (0,1). 

It is easy to show that P^ is identically 1/2 on the 
Po~axis and that it has local minima on the pi-axis at 
Pi = P2 = 1/2 and on the remaining boundary of the 
triangle at Pq = pi — 1/2. There also exists a local 
extrcmum of P^ in the interior of T at the point 
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Fig. 3 (left panel) shows P^ plotted over T for 77 = 
0.6. The interior extremum point {Po,Pi,P2) achieves the 
minimum error probability. For comparison, we consider 
also a signal-only input state of the form 



Ws = ypo|0)s + VPr|l)s + Vp^2\2)s- 



(20) 



For each choice of {po,pi) in T, the error probability 

-prsignal-only , , . ,, . i i - . 

P^ IS computed numerically using (1191), and m 

Fig. 3 (right panel), the difference Pg — Pg is 

plotted on T. The difference is everywhere non-negative, 
being zero on the two boundaries of T other than the 
Pi— axis. 

Lossless Image Sensing: In the lossless case, 1 = with 
probability one, so that the performance of the hypothet- 
ical measurement described in Theorem 1 is attainable 
with any choice of idler states. That performance is de- 



termined by the Gram matrix elements from ( 15 ) 
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where V = {ny->^ . . . , n^^') has praip^. Choosing \(pn)i — 
10)/, the signal-only state 
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FIG. 3. Left; The error probability P^°^ of the NDS state 



of the form of Eq. ( 18 1 for 77 = 0.6. Right: The difference 

Pe ~Pe between the error probabiUties of correspond- 

ing signal-only and NDS states as a function of (po,Pi)- 



with J^P'' = 1, J' = suffices to attain C. In the absence 
of loss, the {/Cm}m=i ^I'S unitary channels. This result 
shows that, among parallel strategies, ancillas do not im- 
prove sensing of M unitary phase images under a signal 
energy constraint. This is unlike the case of minimum er- 
ror probability discrimination of finite-dimensional uni- 
taries in [16], although ancillas are not required for dis- 
criminating two unitaries [ini [U] . The fact that single- 
pass imaging (J'^^^ = 1) suffices is also remarkable as 
there are examples of pairs of unitaries that are better 
(even perfectly) discriminated if multiple shots are al- 
lowed [milH]. 
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APPENDIX: PROOFS OF LEMMAS 1 AND 2 

Lemma 1 (Concavity of C under mixing of ensembles). 
Consider an imaging task with cost function C{m,x). 
For M-ary ensembles {£1} indexed by I, and probability 
distribution {A;}, 
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5] A A >^AiC[£z]. 
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Proof. The LHS of (21) evaluates to 
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The inequality ( 23 ) follows because the inner sum is 



minimized separately for each value oi I in (231 but not 



in (22). D 



Physically, the RHS of (21 1 is the minimum cost 



when one knows the sub-ensemble (indexed by I) prior 



to making the measurement. Then (21) simply asserts 
that the cost is higher in the case where information on 
I is not available. 



Lemma 2 // {£;} are pairwise orthogonal ensem- 



bles on H, we have 



£ = ^; \i£i is then given by 
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Y^XiEi =Y,^iC[£i]- 
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Proof. Let "H; — supp £i for each I. It is clear that it is suf- 
fices to consider POVMs on the subspace %' = 0; T-Li C 
%. For each value of l^ let 11; denote the projection oper- 
ator onto T-Li. Since the {£/} are pairwise orthogonal, we 
have UiHi, = (5;.//n/. Let {E'^'^I^sA' denote the POVM 
(on W) that realizes C [£i]. We define a POVM {E^}:rex 
on %' with elements :- 



s, = ^4')nz, xGA-. 



(25) 



The average cost realized by this POVM on the mixture 
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Here, ( 26 ) follows from 



C{m,x) (26) 

(27) 

(28) 



which in turn follows from the orthogonality of the en- 



sembles {£i}- Equality (27 1 follows from the definition of 
{Ei'^^}. The fact that the POVM {E^} achieves the RHS 



of (21) establishes the claim (24). D 



Physically, the optimum measurement on £ may be re- 
garded as a two-step measurement that learns the value 
of I first and then makes the measurement achieving 
C[£;]. The orthogonality of the {£;} ensures that the 
first measurement can be performed perfectly without 
adding noise to the second. 



